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ABSTRACT 


Calculations  of  turbulent  shear  stress  distributions  are  reported 
for  zero  pressure  gradient  and  adverse  pressure  gradient  supersonic 
turbulent  boundary  layer  flow.  The  calculations  are  accomplished  by 
numerically  integrating  the  equation  for  conservation  of  streamwise 
momentum  using  mean  profile  experimental  data.  The  mixing  length 
distributions  have  also  been  determined  and  indicate  that  the  mixing 
length  distribution  is  significantly  altered  for  the  adverse  pressure 
gradient  flow.  Finite  difference  boundary  layer  computations  using  an 
altered  mixing  length  distribution  show  improved  agreement  with  experi¬ 
mental  measurements  of  skin  friction  for  the  adverse  pressure  gradient 
flow. 
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I .  INTRODUCTION 


Since  no  extensive  direct  measurements  had  been  reported  of  the 
turbulent  shear  stress  distribution  in  the  supersonic  turbulent  boundary 
layer  until  the  recent  experiment  of  Rose1*,  models  of  the  shear  stress 
distribution  have  been  obtained  by  extending  results  found  for  the  more 
thoroughly  measured  incompressible  turbulent  boundary  layer.  This 
extension  is  justified  on  the  basis  of  available  turbulent  fluctuation 
data  that  indicate  little  change  in  the  structure  of  the  turbulent 
fluctuations  for  the  supersonic  (M  =  1. 2-5.0)  turbulent  boundary  layer 
compared  to  that  obtained  in  the  incompressible  turbulent  boundary  layer 
for  zero  pressure  gradient,  adiabatic  flow.  Furthermore,  calculations 
were  reported  by  Maise  and  McDonald2  for  zero  pressure  gradient, 
adiabatic  flow  in  which  the  shear  stress,  eddy  viscosity  and  mixing 
length  distributions  were  determined  by  a  finite  difference  solution  of 
the  boundary  layer  equations  using  a  "law  of  the  wall"  velocity  correla¬ 
tion.  These  calculations  revealed  that  the  mixing  length  distribution 
was  essentially  unchanged  as  the  Mach  number  ranged  from  zero  to  five. 
The  eddy  viscosity  distribution,  however,  was  shown  to  exhibit  a 
sensitivity  to  both  Reynolds  number  and  Mach  number. 

Calculation  procedures  employing  a  shear  stress  distribution 
established  from  incompressible  turbulent  boundary  layer  characteristics 
have,  in  general,  yielded  good  agreement  with  experimental  data  for  the 
supersonic  turbulent  boundary  layer  zero  pressure  gradient  or  mildly 
favorable  pressure  gradient,  adiabatic  flow.  Less  satisfactory  results 
are  obtained  when  computing  flow  in  which  adverse  pressure  gradient  and 
longitudinal  curvature  are  encountered  4»5>6>^  One  obvious  source  for 
error  is  the  inability  of  present  computational  procedures  to  properly 
account  for  the  pressure  gradient  normal  to  the  surface.  Another 
source  of  uncertainty  lies  in  the  behavior  of  the  distribution  of 
turbulent  shear  stress.  Although  turbulent  fluctuation  data  indicate 
similar  behavior  for  incompressible  and  supersonic  zero  pressure 
gradient  turbulent  boundary  layers,  distinct  differences  are  apparent 
in  the  profiles  of  turbulent  fluctuations  between  supersonic  zero 
pressure  gradient  and  supersonic  adverse  pressure  gradient  turbulent 
boundary  layer  flow8’9*  . 

This  report  describes  calculations  of  turbulent  shear  stress, 
mixing  length  and  eddy  viscosity  distributions  using  the  mean  profile 
data  reported  in  reference  9.  These  data  were  obtained  in  a  Mach  3.5 
nozzle  wall  turbulent  boundary  layer  for  zero  pressure  gradient  and  an 
isentropic-ramp-induced  adverse  pressure  gradient  flow  (see  Figure  1). 
The  method  for  calculating  the  shear  stress  distribution  for  the  zero 
pressure  gradient  data  is  similar  to  that  of  Meier  and  Rotta10.  The 
procedure  for  calculating  the  shear  stress  profile  for  the  flow  over 
the  ramp  model  was  developed  from  the  boundary  layer  equations  applica¬ 
ble  to  compressible  flow  over  a  surface  with  longitudinal  curvature. 

*  References  are  listed,  on  page  21, 
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A  similar  procedure  applied  to  experimental  data  for  this  flow  configu¬ 
ration  has  not  been  previously  reported. 

II.  ANALYSIS  OF  ZERO  PRESSURE  GRADIENT  DATA 
A.  Shear  Stress  Profile 


The  equations  for  conservation  of  mass  and  conservation  of  stream- 
wise  momentum  applicable  to  zero  pressure  gradient,  adiabatic  compres¬ 
sible  flow  over  a  flat  plate  can  be  combined  and  integrated  in  the 
direction  normal  to  the  surface  to  yield 
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Assuming  the  flow  to  be  locally  similar  enables  equation  (1)  to  be 
written  in  dimensionless  form  as 
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The  assumption  of  "locally  similar  flow"  (see  references  3,  10  and 
11)  is  not  expected  to  be  fully  satisfied.  However,  in  the  vicinity 
near  the  wall  where  the  assumption  of  similarity  is  poorest,  the  con¬ 
tribution  of  the  convective  terms  is  minimal.  The  equality  (l/6)(d6/dx) 
(l/0)(d9/dx)  also  follows  from  similarity.  Introducing  this  equality 
and  the  relation  (d0/dx  =  cf/2)  into  equation  (2)  yields  after 
rearranging 
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The  integrals  in  equation  (3)  were  evaluated  numerically  using  the 
tabulated  profile  data  of  reference  9.  The  numerical  integration  was 
carried  out  using  the  standard  Fortran  subroutine  AVINT  which  fits  a 
parabola  to  the  data.  An  example  of  the  profile  obtained  is  shown  in 
Figure  2  compared  to  that  obtained  by  Maise  and  McDonald2.  These, 
profiles  are  in  reasonable  agreement  except  near  the  edge  of  the 
boundary  layer.  The  profile  from  reference  2  goes  to  zero  at  y/6  =  1 
with  a  finite  slope  whereas  the  profile  calculated  here  approaches  zero 
asymptotically  beyond  the  edge  of  the  boundary  layer.  This  behavior  is 
consistent  with  the  profile  of  turbulent  fluctuations  as  measured  using 
hot-wire  anemometry  (see  references  8,  9  and  12) . 
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B.  Velocity  Derivative 


Also  needed  for  the  calculation  of  eddy  viscosity  or  mixing  length 
distributions  is  the  distribution  of  the  velocity  derivative,  du/dy. 
This  was  calculated  from  the  experimental  data  using  a  central  differ¬ 
ence  technique.  An  example  of  the  profiles  obtained  is  shown  in  Figure 
3.  No  attempt  was  made  to  smooth  the  data  since  it  was  felt  that  the 
scatter  indicated  by  the  data  would,  in  itself,  be  of  interest  in 
revealing  the  uncertainty  in  the  eddy  viscosity  and  mixing  length 
distributions  obtained.  The  viscous  sublayer  is  extremely  thin 

(y+  =  10  -*■  y  »  .004-inch)  for  the  experimental  data  considered  here. 

The  value  of  d  Cu/u00)/d  (y/S)  corresponding  to  the  measured  wall  shear 
stress  is  approximately  90. 

C.  Eddy  Viscosity  and  Mixing  Length 

The  eddy  viscosity  and  mixing  length  distributions  were  calculated 


using  the  relations 
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where  the  local  value  of  density  was  obtained  from  the  tabulated  profile 
data. 

An  example  of  the  mixing  length  distribution  obtained  is  shown  in 
Figure  4.  The  data  in  the  wall  region  of  the  boundary  layer  agree  very 
well  with  the  accepted  relation  l  =  .4y.  "Hie  magnitude  in  the  region 
of  the  plateau  (about  0.055)  is  considerably  less  than  the  generally 
accepted  value  of  0.089.  The  trend  of  £/6  increasing  for  y/6  >  .6  has 
been  observed  in  other  calculations3  and  will  be  discussed  further  in 
the  next  section. 

An  example  of  the  eddy  viscosity  distribution  obtained  is  shown  in 
Figure  5  compared  to  profiles  from  reference  2.  TTie  peak  value  is 
significantly  less  than  that  in  reference  2.  The  trend  of  the  profile 
for  y/6  >  .6  again  deviates  from  the  trend  of  the  profiles  of  reference 

1. 

D.  Analysis  of  Mixing  Length  Distribution  Anomalous  Behavior 

There  are  two  trends  shown  in  the  mixing  length  distribution  that 
invite  question:  (1)  the  low  value  of  the  plateau  region,  and  (2)  the 
behavior  for  y/6  >  0.6.  The  cause  for  the  seemingly  anomalous  behavior 
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was  sought  by  comparing  the  results  obtained  when  the  shear  stress 
distribution  calculated  here  was  replaced  by  that  of  reference  2,  and 
likewise,  for  the  distribution  of  the  velocity  derivative,  du/dy. 

The  mixing  length  distribution  obtained  using  the  shear  stress 
profile  of  reference  2  was  virtually  identical  to  that  in  Figure  4, 

The  next  step  taken  was  to  alternate  profiles  of  the  velocity  derivative. 
The  velocity  profile  used  in  reference  2  was  a  compressible  "law  of  the 
wall"  correlation.  Equation  11  of  reference  2,  which  was  obtained  by 
inverting  the  velocity  defect  correlation  is  repeated  below. 


where 


2-5  loge  ($] 


+  1.25 


l1  +  cos  Hi) 


+  arcs in 


[Cy-D/2]  M l 
l  +  [Cy-D/2]  m l 


(6) 


This  expression  was  differentiated  to  yield  the  velocity  derivative 
in  the  y  direction.  The  constant  1.25  in  equation  (6)  is  the  ratio, 
n/K,  where  n  is  Coles*  wake  parameter  and  K  is  the  constant  in  the 
relation  £  =  Ky.  Using  the  calculated  velocity  derivative,  the  density 
profile  from  the  tabulated  data,  and  both  the  shear  stress  profile 
calculated  here  and  that  from  reference  2,  the  value  of  £/6  in  the 
plateau  region  was  about  0.11. 

Since  a  value  of  n/K  of  2.25  was  found  to  correlate  the  experimental 
data  of  reference  9,  this  value  was  substituted  into  equation  (6).  The 
mixing  length  distribution  now  obtained  is  shown  in  Figure  6  for  the  two 
shear  stress  distributions.  The  value  of  £/6  in  the  plateau  region  is 
about  0.070.  The  trend  for  £/6  to  increase  for  y/6  >  0.6  is  obviously 
the  effect  of  the  calculated  shear  s  tress  profile  approaching  zero 
asymptotically  beyond  y/6  =  1.0. 

In  comparing  the  velocity  distribution  calculated  using  equation 
(6)  for  n/K  =  2.25  with  the  experimental  data,  a  significant  disagree¬ 
ment  is  noted  in  the  region  for  y/6  <  0.4.  This  is  shown  in  Figure  7. 

The  effect  of  this  discrepancy  on  the  velocity  derivative  is  shown  in 
Figure  8.  The  region  of  significant  disagreement  lies  between  y/6  =  0.15 
and  y/6  =  0.6.  This  is  the  region  in  which  the  plateau  forms  as  noted 
in  Figure  4.  Although  the  data  correlate  in  velocity  defect  coordinates, 
this  type  of  correlation  is  not  sufficiently  accurate  to  extract  a  valid 
profile  for  the  velocity  derivative. 
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The  conclusions  reached  regarding  the  mixing  length  distribution 
are:  (1)  the  profile  of  Figure  4  is  a  valid  representation  of  the 
effective  mixing  length  for  the  experimental  data  considered;  and  (2)  a 
"law  of  the  wall"  velocity  correlation  is  not  a  sufficiently  accurate 
representation  of  the  velocity  profile  to  enable  a  valid  velocity 
derivative  to  be  obtained. 


III.  ANALYSIS  OF  ADVERSE  PRESSURE  GRADIENT  DATA 
A.  Shear  Stress  Distribution 


The  equations  for  conservation  of  mass  and  streamwise  momentum  as 
applicable  to  two-dimensional,  compressible  turbulent  boundary  layer 
flow  over  a  surface  with  longitudinal  curvature  may  be  expressed  as 
follows. 
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These  relations  result  from  order  of  magnitude  considerations 
applied  to  the  more  general  relations  derived  by  Tetervin13. 

Equations  (7)  and  (8)  may  be  integrated  in  the  y  direction  normal 
to  the  local  surface  and  combined  to  yield  the  following  relation  for 
the  shear  stress  distribution. 
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where 


1 

1+icy 


(3  = 


and  t  =  y  -  pi/  v*  . 


This  relation  has  been  evaluated  directly  using  measured  values  of  wall 
shear  stress  and  the  tabulated  profile  data  of  reference  9.  Values  of 
the  partial  derivatives  have  been  determined  along  lines  of  constant 
mass  flux  using  a  least  squares  technique  which  fits  a  parabola  to  the 
data.  The  lines  of  constant  mass  flux  are  shown  in  Figure  9  along  with 
the  change  in  boundary  layer  thickness.  It  is  interesting  to  note  that 
although  the  boundary  layer  becomes  less  thick  at  downstream  stations 
mass  continues  to  be  entrained  within  the  boundary  layer.  An  example 
of  the  distribution  of  the  partial  derivative  3(pu)/3x  is  shown  in 
Figure  10.  In  evaluating  the  shear  stress  distribution,  equation  (9) 
was  integrated  numerically  as  discussed  previously. 


Initial  efforts  to  calculate  the  shear  stress  distribution  resulted 

in  profiles  with  a  large  (about  6  t  )  negative  value  in  the  vicinity  of 

•  w 

the  boundary  layer  outer  edge.  After  examining  the  accuracy  of  the 
profiles  of  the  partial  derivatives,  it  was  concluded  that  these  partial 
derivatives  could  be  in  error  by  as  much  as  ±  30  percent.  In  order  to 
obtain  physically  meaningful  profiles  for  the  shear  stress  distribution, 
the  partial  derivative  of  mass  flux,  3(pu)/3x,  was  corrected  as  indica¬ 
ted  below.  The  corrected  profile  is  also  shown  in  Figure  10. 
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A  computation  iteration  loop  was  used  to  determine  the  value  of  the 
constant  A  that  resulted  in  the  shear  stress  profile  remaining  greater 
than  zero  throughout  the  boundary  layer.  The  value  of  A  varied  from 
station  to  station  and  ranged  between  0.7750  at  the  first  station  and 
0.6875  at  the  last  station  on  the  ramp  model.  The  shear  stress  pro¬ 
files  obtained  are  shown  in  Figure  11.  These  profiles  are  consistent 
with  the  boundary  condition  at  the  wall,  dx/dy  =  dp/dx,  and  with  the 
trend  of  the  fluctuation  data  shown  in  Figure  12.  The  turbulent 
fluctuation  data  were  obtained  using  constant  temperature  hot  wire 
anemometry.  These  data  are  plotted  as  an  arbitrary  dimensionless  number 
determined  by  dividing  the  measured  signal  by  0.005  volt.  Returning  to 
Figure  11,  the  trend  of  the  shear  stress  profile  to  increase  beyond  the 
point  of  minimum  shear  stress  is  inconsistent  with  the  requirement  that 
x  approach  zero  in  the  vicinity  of  the  edge  of  the  boundary  layer.  This 
trend  is  believed  to  be  a  result  of  the  uncertainty  in  determining  the 
streamwise  partial  derivatives. 
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A  check  was  made  to  assess  the  effect  on  the  shear  stress  profiles 
of  adjusting  the  partial  derivative  3(pu2)/3x  and  comparing  the  results 
with  those  obtained  by  correcting  3(pu)/9x.  The  profiles  were  in  close 
agreement  both  in  trend  and  magnitude. 

B.  Velocity  Derivative 


The  velocity  derivative,  du/d y,  was  calculated  directly  from  the 
tabulated  profile  data  in  the  same  manner  as  discussed  previously.  An 
example  of  the  velocity  derivative  for  the  flow  over  the  ramp  model  is 
given  in  Figure  13.  This  profile  is  significantly  unlike  that  of  Figure 
3  in  the  wake  portion  of  the  boundary  layer.  The  greater  value  of  du/dy 
in  this  portion  of  the  boundary  layer  is  consistent  with  the  trend  of 
the  fluctuation  data  since  (I?  2)  can  be  considered  to  be  proportional  to 
(du/dy)2. 

C.  Mixing  Length  and  Eddy  Viscosity  Profiles 

The  distributions  of  mixing  length  and  eddy  viscosity  were  calcu¬ 
lated  using  equations  (4)  and  (5)  as  discussed  previously. 

An  example  of  the  mixing  length  profiles  obtained  is  shown  in 
Figure  14.  These. profiles  are  considerably  distorted  compared  to  that 
obtained  in  the  region  of  zero  pressure  gradient.  The  first  character¬ 
istic  of  interest  is  the  change  in  value  of  K  in  the  relation  i  =  Ky. 

A  value  for  K  of  0.65  appears  to  fit  the  data  points  in  the  vicinity 
of  the  wall.  A  second  characteristic  of  interest  is  the  increase  in 
peak  value  at  succeeding  downstream  stations.  The  behavior  of  the 
mixing  length  profiles  for  y/<5  >  0.5  is  suspect  due  to  the  uncertainty 
associated  with  the  calculated  shear  stress  profiles. 

The  distribution  of  eddy  viscosity  is  shown  in  Figure  15.  These 
profiles  are  similar  to  the  profiles  obtained  for  zero  pressure  gradient 
flow  with  the  following  exceptions:  (1)  the  peak  value  changes  at 
different  streamwise  positions;  and  (2)  the  magnitude  of  the  peak  value 
is  from  2  to  4.5  times  the  value  obtained  for  the  zero  pressure 
gradient  flow. 


IV.  NUMERICAL  CALCULATION  USING  MODIFIED  MIXING  LENGTH  DISTRIBUTIONS 

The  boundary  layer  properties  have  been  calculated  using  the 
computer  program  of  Hixon,  Beckwith  and  Bushnell15.  This  program  is  an 
implicit  finite-difference  procedure  that  calculates  turbulent  flow 
using  an  eddy  viscosity  model  based  on  a  mixing  length  distribution. 

The  mixing  length  distribution  is  modified  for  the  effect  of  wall 
damping  in  a  very  small  region  close  to  the  wall  by  using  Van  Driest* s 
exponential  damping  function.  Calculations  have  been  run  using  the 
conventional  mixing  length  distribution  as  well  as  those  calculated 
here. 
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The  conventional  mixing  length  distribution  used  was: 


Jt/6  =  0.4  y/6 


Jl/6  =  0.04  +  ^  (y/6)  ~  0,1  j  (. 


Jl/6  =  0.10 


for  y/6  <0.1 

06)  for  0.1  <  y/6  <  0.3  (12) 

for  y/6  >  0.3 


Other  mixing  length  distributions  used  were: 

(I)  zero  pressure  gradient 

Jt/6  =  0.4  y/6  for  y/6  <0.1 


A/6  = 

0.04  ♦[(y/6,.2 

°*1J  (0.015) 

for 

0. 1  <  y/6  <0.3 

Jt/6  =  .055 

for 

0.3  <  y/6  <  0.65 

Jl/6  = 

.1857  (y/6  -  .65)  +  .055  for 

0.65  <  y/6  <  1 

Jt/6  =  .12 

for 

y/6  >  1  . 

(2)  adverse  pressure  gradient 

(a) 

Jt/6  =  0.65  y/6 

for 

y/6  <0.1 

Jt/6  -  0.065  + 

fcy/o  -  Q-O 

(.035) 

for  0.1  <  y/6  <  0.3 

Jt/6  =  .08 

for 

y/6  >  0.3 

(b) 

Jt/6  =  0.5  y/6 

for 

y/6  <0.1 

A/6  =  0.05  + 

r  (y/6)  -  0.1  1 

L  -2  J 

(0.05) 

for  0.1  <  y/6  <  0.3 

00 

O 

• 

II 

to 

\ 

o* 

for 

y/6  >  0.3 

(13) 


(14) 
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The  calculations  were  started  using  measured  profiles  of  velocity 
and  total  temperature.  The  program  then  calculates  the  density  profile 
assuming  constant  static  pressure  across  the  boundary  layer.  The 
effective  length  or  starting  point  for  the  boundary  layer  development 
has  been  taken  as  the  nozzle  throat  for  the  calculations  shown  here. 

Velocity  profiles  for  the  zero  pressure  gradient  flow  calculated 
using  the  mixing  length  distributions  given  as  equations  (12)  and  (13) 
are  shown  in  Figure  16  compared  to  the  experimental  data.  The  calcula¬ 
tions  using  the  conventional  mixing  length  distribution  predict  a 
velocity  profile  that  is  too  full  compared  to  the  measured  profile; 
however,  the  calculations  using  the  mixing  length  distribution  determined 
here  (equation  13)  show  very  good  agreement  with  the  experimental  data. 

It  cannot  be  concluded  that  these  results  are  true  for  the  compressible 
turbulent  boundary  layer  in  general.  However,  the  response  of  the 
numerical  program  to  the  calculated  mixing  length  distribution  does 
indicate  that  equation  (13)  is  a  valid  representation  of  the  effective 
mixing  length  distribution  for  the  experimental  data  considered. 

Calculated  results  for  skin  friction  coefficient  are  shown  in 
Figure  17  compared  to  experimental  measurements  of  skin  friction  coef¬ 
ficient.  The  conventional  mixing  length  distribution  yields  good 
agreement  with  the  experimental  data  for  zero  pressure  gradient  flow. 

The  agreement  for  the  calculated  mixing  length  distribution  is  also 
within  the  accuracy  of  the  experimental  data.  As  the  calculation 
proceeds  downstream,  the  skin  friction  coefficient  begins  to  decrease 
when  the  flow  encounters  the  adverse  pressure  gradient.  Values  for 
skin  friction  coefficient  calculated  using  the  conventional  mixing 
length  distribution  are  15-20%  low  compared  to  the  experimentally 
measured  values  for  the  adverse  pressure  gradient  flow. 

Additional  calculations  for  the  adverse  pressure  gradient  flow  are 
shown  in  Figure  18.  These  calculations  were  started  at  the  first  test 
station  on  the  ramp  model  using  measured  profiles  of  velocity  and  total 
temperature.  Again,  the  calculated  values  of  skin  friction  coefficient 
are  about  20%  low  compared  to  the  experimental  measurements  when  the 
conventional  mixing  length  distribution  is  used.  The  calculations  using 
a  value  for  K  of  0.65  yield  values  of  skin  friction  coefficient  that  are 
about  40%  too  great.  However,  the  calculations  run  using  a  value  for  K 
of  0.50  do  show  promise.  These  results  agree  within  ±  5%  of  the 
measured  values. 


V.  CONCLUSIONS 

Calculations  have  been  made  of  the  turbulent  shear  stress  distri¬ 
bution  for  compressible  zero  pressure  gradient  and  adverse  pressure 
gradient  boundary  layer  flow.  The  effective  mixing  length  and  eddy 
viscosity  distributions  have  been  determined  using  these  calculated 
distributions  of  turbulent  shear  stress. 
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In  order  to  obtain  physically  meaningful  distributions  of  turbulent 
shear  stress  for  the  adverse  pressure  gradient  flow,  it  was  necessary  to 

adjust  the  profile  for  the  partial  derivative,  (pu) ,  that  was  deter¬ 
mined  directly  from  the  tabulated  profile  data.  It  is  recognized  that 
the  resulting  shear  stress  profiles  invite  some  uncertainty  due  to  this 
procedure.  However,  it  is  felt  that  the  calculations  reported  are 
physically  meaningful  and  that  the  results  indicate  trends  which  could 
lead  to  improved  prediction  of  compressible  turbulent  boundary  layer 
development  in  an  adverse  pressure  gradient. 

The  conclusions  reached  on  the  results  of  this  investigation  are 
summarized  below. 

(1)  The  mixing  length  distribution  shown  in  Figure  4  is  a  valid 
representation  of  the  effective  mixing  length  distribution  for  the  zero 
pressure  gradient  data  considered. 

(2)  The  value  of  the  constant  K  in  the  relation  £  =  Ky  changed 
from  0.4  for  zero  pressure  gradient  flow  to  0.65  for  the  adverse 
pressure  gradient  flow. 

(3)  The  maximum  values  of  t/t^  increased  as  the  flow  proceeded  in 

the  streamwise  direction  in  the  region  of  adverse  pressure  gradient. 

This  increase  in  (t/t  )  resulted  in  the  value  of  (l/S)  likewise 
v  w  max  'max 

increasing  at  downstream  stations  in  the  adverse  pressure  gradient. 

(4)  Both  the  conventional  and  the  calculated  mixing  length 
distributions  yielded  good  agreement  with  experimental  measurements  of 
wall  shear  stress  for  zero  pressure  gradient  flow  when  used  in  a  finite 
difference  boundary  layer  computer  program. 

(5)  The  best  agreement  of  the  numerical  calculations  with  experi¬ 
mental  measurements  of  wall  shear  stress  for  the  adverse  pressure 
gradient  flow  was  obtained  by  substituting  K  =  0.5  in  the  mixing  length 
distribution. 

Since  the  computer  program  used  for  these  numerical  calculations 
is  restricted  to  a  constant  static  pressure  through  the  boundary  layer, 
it  would  be  of  interest  to  perform  calculations  of  the  flow  over  the 
ramp  model  using  a  program  that  accounts  for  the  pressure  gradient 
normal  to  the  surface. 
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Figure  1.  Tunnel  Wall  Configuration  and  Test  Station  Identification 


Figure  2.  Turbulent  Shear  Stress  Distribution 


24 


t - r 


Moo=3.53 
dp /dx  =  0 


T 


.6 

.4 


.2 

0 


J _ i  f _ | _ L 

1.0  2.0 

d(u/Uoo) 

d(y/8) 


Figure  3.  Velocity  Derivative  Profile,  dp/dx  =  0 
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Experimental  Velocity  Profile  Compared 
to  the  Law  of  the  Wall  Correlation 
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Figure  8.  Experimental  Velocity  Derivative  Compared 
to  the  Result  Obtained  Using  the  Law  of 
the  Wall  Correlation 
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Figure  9.  Lines  of  Constant  Mass  Flux  for  Flow  Over  the  Ramp  Model 


Figure  10.  Distribution  of  the  Streamwise  Derivative 
of  Mass  Flux 
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Figure  12.  Turbulent  Fluctuation  Data  Obtained  Using 
Constant  Temperature  Hot  Wire  Anemometry 
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Figure  13.  Velocity  Derivative  Profile,  dp/dx  >  0 
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Figure  17.  Numerical  Calculations  of  Skin  Friction 

Coefficient  Compared  to  Experimental  Measurements 


STREAMWISE  DISTANCE , INCHES 
117  118  119  120  121  122 

TEST  STATION  DESIGNATION 


Figure  18.  Comparison  of  Skin  Friction  Calculations  for  Adverse 

Pressure  Gradient  Flow  Using  Unconventional ^Values  for  K 
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Calculations  of  turbulent  shear  stress  distributions  are  reported  for  zero 
pressure  gradient  and  adverse  pressure  gradient  supersonic  turbulent  boundary 
layer  flow.  The  calculations  are  accomplished  by  numerically  integrating  the 
equation  for  conservation  of  streamwise  momentum  using  mean  profile  experimental 
data.  The  mixing  length  distributions  have  also  been  determined  and  indicate 
that  the  mixing  length  distribution  is  significantly  altered  for  the  adverse 
pressure  gradient  flow*  Finite  difference  boundary  layer  computations  using  an 
altered  mixing  length  distribution  show  improved  agreement  with  experimental 
measurements  of  skin  friction  for  the  adverse  pressure  gradient  flow. 
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turbulent  shear  stress  distribution 
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supersonic  turbulent  boundary  layer 
zero  pressure  gradient 
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